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Abstract. We study theoretically light propagations at the zigzag edge of a 
honeycomb photonic crystal consisting of dielectric rods in air, analogous to graphene. 
Within the photonic band gap of the honeycomb photonic crystal, a unimodal edge 
state may exist with a sharp confinement of optical fields. Its dispersion can be tuned 
simply by adjusting the radius of the edge rods. For the edge rods with a graded 
variation in radius along the edge direction, we show numerically that light beams of 
different frequencies can be trapped sharply in different spatial locations, rendering 
wideband trapping of light. 



PACS numbers: 42.25.Bs, 42.25.Gy, 42.70. Qs 



Submitted to: J. Phys.: Condens. Matter 



Wideband trapping of light by edge states in honeycomb photonic crystals 



2 



Graphene, a flat monolayer of graphite with carbon atoms arranged in a honeycomb 
lattice, has many interesting physical properties with great promise in nanoelectronics 
owing to its unique electronic properties, e.g., conical dispersion (Dirac cones) in the 
electronic band structure [1]. Inspired by graphene, its photonic analog, namely, two- 
dimensional photonic crystals [2H1] that possess Dirac cones in photonic band structures, 
have drawn considerable attention in recent years [514T2] . The existence of Dirac cones 
may lead to many interesting optical phenomena such as conical diffraction and gap 
sohtons [5], pseudo diffusive transport of light |6l|TT], zitterbewegung of photons [8] and 
the realization of zero-refractive-index materials [T2] . 

Wideband trapping of light, also called "Rainbow" trapping, is an optical 
phenomenon that light of different frequencies can be slowed down and stopped 
at different spatial positions in specially designed waveguides [13]. It has many 
potential applications ranging from optical information storage to enhanced light-matter 
interactions [H]. Several structures for reahzing wideband trapping have been proposed, 
including metamaterial waveguides [131 [15], plasmonic waveguides [16], metal-coated 
optical nano waveguides [17] and plasmonic gratings [l8l[T9] . However, these structures 
may suffer inevitably a serious problem of optical loss due to constituting metals, which 
are rather lossy especially in the visible regime. In contrast, dielectric photonic crystal 
waveguides can also slow light with a very small loss [20l|21] , which could be exploited in 
light trapping [22]. But they are usually multimodal [21], which may lead to undesirable 
mode coupling. 

By breaking the inversion symmetry of photonic crystals that possess Dirac cones, a 
photonic band gap (PBG) around Dirac points may open up [23]. With the presence of 
edges, edge states within the PBG may appear with a sharp confinement of optical fields 
at the edges [24|. With further breaking the time- reversal symmetry, confined edge states 
become even unidirectional, showing robustness against scattering or disorder [7t [25] - [28] . 
In this paper, we show theoretically that edge states in honeycomb photonic crystals 
(HPCs) can be exploited in realizing wideband trapping of light. Since no metal is 
introduced, the optical loss should be much lower than those in metamaterial and 
plasmonic waveguides. Compared with photonic crystal waveguides, our system is 
unimodal and more compact as well. Although exemplified by HPCs with rods being 
infinitely long, our idea can be applied to two-dimensional photonic crystals with finite 
rod lengths or membrane photonic crystals with drilled holes [201422] . 

HPCs under study consist of dielectric rods which are arranged in a honeycomb 
lattice in air, analogous to graphene, as schematically shown in figure lU^a). The inter- 
distance between adjacent rods is denoted by d. In an infinite HPC, there are two 
different sets of rods whose radii are denoted respectively by vi and r2. The radius of 
edge rods, which are changeable, is denoted by r^. The lattice constant of an HPC is 
thus a = \/3d. Without loss of generahty, the dielectric constant of all rods is taken to 
be 11.4, a typical value for silicon at telecommunication wavelengths. In our following 
discussions, we only consider the transverse magnetic polarization, i.e., the case that 
the electric field is parallel to the rods. 
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Figure 1. (a) Top view of an HPC witli a zigzag edge. Light and dark gray solid circles 
represent different sets of rods in a honeycomb lattice. Edge rods are marked by yellow 
solid circles. Doted lines represent the unit cell of the HPC. (b) First Brillouin zone 
of the honeycomb lattice with F, M and K representing points of high symmetry, (c) 
Photonic band structures of HPCs. Red lines are those for ri — r2 — OAd, while blue 
lines for ri = OAd and r2 = 0.2d. The colored area stands for the PEG of the latter 
HPC. (d) PEG (colored area) map as a function of r2 for fixed ri — OAd. Frequency 
is in units of a/ A, where a is the lattice constant and A is the wavelength in vacuum. 

In figure [U^c), photonic band structures of infinite HPCs, calculated by a plane- 
wave expansion method [1], are shown. For an HPC possessing the inversion symmetry, 
i.e., Ti = r2, Dirac points exist where the first and second photonic bands cross each 
other at the K point with a reduced frequency of a/A = 0.253. With broken inversion 
symmetry, i.e., vi ^ r2, the degeneracy at the K point will be lifted and a PEG will open 
up [23]. The lower (upper) edge of the PEG is defined by the first (second) photonic 
band at the K (M) point. 

To attain an HPC with an optimal or desired PEG, one can tune the radius of the 
rods or adjust the dielectric constant of the rods. In figure [T]^d), the PEG map as a 
function of r2 for fixed ri = OAd is shown. For ri = r2 = OAd, there is no PEG due to 
the existence of the Dirac points. For ri 7^ r2, a PEG always exists. For ri = OAd and 
r2 = 0.2d, the PEG ranges in reduced frequency a/A from 0.26 to 0.343. For ri = OAd 
and r2 = 0, the relative PEG width, Au/u^, can be as high as 43%, where Au is the 
PEG width and cUm is the mid-gap frequency. 

Figure M^a) shows the dispersion of edge states of an HPC with a zigzag edge, 
calculated by the plane- wave expansion method [4j . In our calculations, a supercell with 
a size of a along the edge direction and AGd along the direction perpendicular to the edge 
direction is used. Periodical boundary conditions are imposed in both directions. The 
supercell contains a finite HPC with a thickness of 35d along the direction perpendicular 
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Figure 2. (a) Calculated dispersion of edge states for an HPC with a zigzag edge as 
shown in figure HJa) , where fc| | is the wavevector along the edge direction in units of 
n/a. The radius parameters for the HPC are ri = OAd and r2 = 0.2d. Gray areas 
represent projected photonic bands. Gray, blue, and dark gray lines represent the edge 
states for the edge rods with different radius. Dotted straight lines stand for the light 
line in vacuum. Inset shows the cutoff frequency (at fc|| = 7r/a) as a function of r^. 
Two doted lines defines the range of the PEG. (b)-(d) Distributions of the electric 
field intensity \E\'^ for the edge state with — 0.28d at (b) fc|| = O-Gn/a, (c) O.Sir/a, 
and (d) 7r/a, marked by dots in (a). Doted circles stand for the rods. 



to the edge direction, corresponding to 48 rows of rods including the edge rods. Since the 
separation of the two zigzag edges are large enough, there is no coupling between the two 
edges. The resulting dispersion can thus be considered as that for a single zigzag edge 
of a semi-infinite HPC. Projected photonic bands are obtained by projecting the bulk 
photonic band structure over the wavevector along the edge direction, k^^. In addition 
to photonic bands, for a given radius of the edge rods r^, an edge state may appear 
between the first and second photonic bands. Interestingly, this edge state is unimodal, 
whose dispersion depends on r^. 

Even for different values of r^, the dispersion of the edge state has a similar 
behaviour: its frequency increases with increasing k\\ up to the Brillouin zone boundary 
k\\ = n/a. At fell = vr/a, it takes a zero slope. In other words, the group velocity 
decreases with increasing k\\ and approaches zero at the Brillouin zone boundary. This 
interesting feature could be exploited to slow down and even stop light |14l|20]. The 
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frequency of an edge state at A;|| = 7r/a is dubbed cutoff frequency since the frequencies 
of the edge state at other wavevectors are always smaller than the cutoff frequency. 
For different values of r^, the corresponding cutoff frequency is different, showing a 
monotonic decrease with increasing r^, as shown in the inset of figure [2]^a). 

Note that edge states lie below the light line, implying that they are evanescent 
waves. In Figs. [2](b)-(d), the field distributions for the edge state with = 0.28d at 
different k\\ are shown. Obviously, optical fields are confined to the edge rods. For the 
edge mode at A;|| = 0.67r/a, it is near the first photonic band, leading to some bulk 
excitations in the vicinity of the edge, which can be clearly seen from figure EJ^b). In 
contrast, for the edge modes away from the photonic bands there are nearly no bulk 
excitations. As shown by us |25, these edge modes are of subwavelength, similar to 
surface plasmons bounded at metallic surfaces [29]. 

As aforementioned, an edge state at its cutoff frequency has a zero group velocity 
and the cutoff frequency decreases monotonically with increasing r^. These features can 
render wideband trapping of light by taking advantage of edge states, provided that the 
edge rods are designed to have a spatially graded variation in radius. If the grade is small 
enough, the graded edge can be considered approximately as a series of small edges, in 
each of which the edge rods possess a constant radius. As a result, each location seems 
to have a different cutoff frequency. For a wave propagating along such a graded edge, 
its group velocity will be gradually reduced along the edge direction and eventually 
approach zero in a location where the corresponding cutoff frequency coincide with the 
frequency of the wave. In other words, the wave will stop around this location. 

Figure |3] shows light trapping for different frequencies at the graded zigzag edge of 
an HPC, simulated by a finite-difference time-domain (FDTD) method [30j. The radii 
of the edge rods are assumed to vary linearly with the position of the edge rods, namely, 
r3 = 0.3d -|- 5 X 10~^x, where x is the position along the edge direction. In our FDTD 
simulations, a continuous Gaussian beam with an intensity of 1 and a width of 2d is 
launched at x = —46a. Perfectly matched layer absorbing boundaries [31] are used. 

Obviously, from figure |3] light beams of different frequency can be trapped spatially 
in different locations. The scenario is as follows. For a right-propagating Gaussian 
beam, it will propagate rightwards after launching. The beam will be slowed down 
towards right since its group velocity decreases rightwards due to the reduction of the 
radius of the edge rods. It will eventually be stopped around the location where the 
corresponding cutoff frequency of the location coincides with the beam frequency. By 
inspecting the beam with a reduced frequency of 0.286, for example, the maximum 
intensity of the standstill beam occurs around the position x ~ —4a. The radius of the 
edge rods around this location is about 0.297(i, very close to but a bit smaller than the 
radius of 0.3d that offers a cutoff frequency of 0.286 [see the inset of figure [2]( a)]. For 
the standstill beam, there are some small additional peaks on the left side of the main 
peak. This is because the right-propagating wave may get reflected more or less since 
the grade of the edge rods is not small enough. The reflected and the right-propagating 
waves may interfere with each other, leading to these small additional peaks, which 
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Figure 3. FDTD simulations of tlie electric field intensity \E\'^ for an HPC with 
a graded zigzag edge at different reduced frequencies (a) 0.276, (b) 0.286, and (c) 
0.297. The radius parameters for the HPC are ri = OAd and r2 — 0.2d. A right- 
propagating Gaussian beam of different frequencies is launched at x = —46a. Circles 
denote dielectric rods. The bottom (top) horizontal axis indicates the position (radius 
of edge rods) along the edge direction. 



could be suppressed in principle by using a finer grade. 

In summary, we studied theoretically the edge states at the zigzag edge of HPCs 
consisting of dielectric rods in air for the transverse magnetic polarization. With broken 
inversion symmetry, a PBG exists between the first and second photonic bands. A uni- 
modal edge state may appear within the PBG, whose dispersion and cutoff frequency 
can be tuned by changing the radius of the edge rods. For a zigzag edge with a gradually 
graded radius in edge rods, light beams of different frequencies can be slowed down and 
trapped spatially in different locations. Our results show that edge states in HPCs could 
be exploited in optical storage and enhanced light-matter interactions. 

This work was supported by the 973 Program (Grant nos. 2011CB922004, 
2012CB921604 and 2013CB632701) and the NSFC. 
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